D-D MIXING AND NEW PHYSICS: 
GENERAL CONSIDERATIONS AND CONSTRAINTS ON THE MSSM 
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Combining the recent experimental evidence of D-D mixing, we extract model-independent in- 
formation on the mixing amplitude and on its CP-violating phase. Using this information, we 
present new constraints on the flavour structure of up-type squark mass matrices in supersymmetric 
extensions of the Standard Model. 



The study of meson oscillations represents one of the 
most powerful probes of New Physics (NP) currently 
available. The K and Bd systems are very well studied 
experimentally and all the measurements performed up 
to now are compatible with the Standard Model (SM) ex- 
pectation, although there is still room for NP which could 
be revealed with improved theoretical tools and experi- 
mental facilities hopefully available in the future [l|, l2| . 

As far as the B s is concerned, the experimental evi- 
dence of oscillation was found only recently at the TeVa- 
tron [3j. While the oscillation frequency is already very 
well known, information on the phase of the mixing am- 
plitude is still quite vague, leaving ample room for exper- 
imental improvements expected from hadronic colliders. 

All this experimental information allows to put model- 
independent constraints on NP contributions to the mix- 
ing amplitudes involving down-type quarks These 
constraints already induce highly non-trivial bounds on 
the flavour structure of many extensions of the SM. 
In particular, considering the Minimal Supersymmetric 
Standard Model (MSSM), the flavour properties of the 
down-type squark mass matrices have already been thor- 
oughly analyzed [3]. 

On the other hand, up to now no evidence was found 
of oscillations of mesons involving up- type quarks. Cor- 
respondingly, the off-diagonal entries of up-type squark 
mass matrices were only weakly bounded [!, @- It is 
remarkable that one of the proposed mechanisms to ex- 
plain the flavour structure of the MSSM and to suppress 
the unwanted SUSY contributions to Flavour-Changing 
Neutral Current (FCNC) processes, namely alignment of 
quark and squark mass matrices [7|, naturally produces 
sizable effects in the up-type sector. In the absence of 
stringent experimental information, these models were 
not tightly constrained JfJ . 

Very recently, BaBar Q and Belle independently 
reported evidence for D-D mixing. In this letter we 
use this information, combined with previous constraints 
on D mixing llj, [PJ, Ilia, LL4| > to put model-independent 
bounds on the mixing amplitude and to constrain the 



relevant entries of the up- type squark mass matrices. 
To fulfil this task we use the mass-insertion approxima- 
tion. Treating off-diagonal sfermion mass terms as inter- 
actions, we perform a perturbative expansion of FCNC 
amplitudes in terms of mass insertions. The lowest non- 
vanishing order of this expansion gives an excellent ap- 
proximation to the full result, given the tight experimen- 
tal constraints on flavour changing mass insertions. It 
is most convenient to work in the super-CKM basis, in 
which all gauge interactions carry the same flavour de- 
pendence as SM ones. In this basis, we define the mass 
insertions {S'i 2 ) AB as the off-diagonal mass terms con- 
necting up-type squarks of flavour u and c and helicity 
A and B, divided by the average squark mass. 

Let us first discuss the recent experimental novelties. 
BaBar studied D° -> K + -k~ and D° -> R-tt+ decays 
as a function of the proper time of the D mesons. As- 
suming no CP violation in mixing, which is safe in the 
SM, this analysis allows to measure the parameters x' 2 
and y' , defined in terms of the mixing parameters x and 
y through the relations: 

x = xcosSKTr+ysinSiiTr, y = — x sixiSkw +y cos Sk w , 

where Skk is the relative strong phase between the 
Cabibbo-suppressed D° — ► K + n~ decays and the 
Cabibbo-favoured D° — > K~tt + ones^ This phase has 
been recently measured by CLEO-c ll]. From a fit to 
D° and D° decays, BaBar is able to exclude the point 
x = y' = (which corresponds to the no-mixing sce- 
nario) with a 3.9er significance (including systematic ef- 
fects). In addition, the BaBar collaboration fitted sepa- 
rately the parameters x'± and y' ± of D° — > K ± 7r T decays 
allowing for CP violation. 

Belle directly determines x and y, studying the D a — > 
7r~ Dalitz plot. In this way, one can separately 
measure the mixing parameters and the strong phase. 
Even though this analysis is not precise enough to claim 
the observation of D-D mixing, it allows to disentan- 
gle mixing parameters from the strong phase 5ktt, when 
D° — » Kg7r^~ tt~ and D° — > Ktt results are combined. 
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Belle also found evidence of D-D mixing, observing 
a deviation from zero (at 3.2,7 including the systematic 

error) of yep = ^~pPZ7fc~] — 1 an d in addition mea- 
sured the CP asymmetry ^ r = (r(D -» KK) - T{D — » 
KK))/(r(D -> ifif) + r(L> -» i^i^)). 

We assume that CP violation can occur in mixing but 
not in decay amplitudes, since the latter are dominated 
by SM tree-level contributions. Therefore, we assume 
that Ti2 is real. Our aim is to determine the parame- 
ters |Mi2|e _i * 12 and Ti2 from the available experimental 
data. One can write the following relations 15]: 
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where S = \p\ 2 — \q\ 2 and <f> is the phase of the mixing 
parameter q/p. We fit for |Mi 2 |, |r i2 1 and $i 2 using the 
experimental inputs listed in Table HI taking into account 
the correlations between y'± and x± in the BaBar results. 
Notice that all observables can be written in terms of 
|M 12 |, |r 12 | and $i 2 . 

The results of the simultaneous fit are quoted in Tab.HTl 
and shown in Fig. [TJ Our results are slightly different 
from the HFAG May 2007 averages because they are ob- 
tained allowing for CP violation, while the HFAG results 
assume no CP violation. In Fig. [21 the two-dimensional 
constraints on the y vs x, 4> vs \q/p\ and <&i 2 vs |Mi 2 | are 
given. We notice that, since the measured value of yep 
is large, the phase <j> is constrained to be close to zero. 
However, due to the large value of T 12 , the constraint 
on $ 12 is less stringent. Some of the results collected 
in Tab. [TT] can be compared with the existing literature. 
Concerning the upper bound on |Mi 2 |, we find an im- 
provement of almost an order of magnitude with respect 



to the analysis of ref. 15j , while for x the improvement 
with respect to ref. [l7| is about a factor of three. 

The calculation of |Mi 2 | is plagued by long-distance 
contributions [18(. To take them into account, we pro- 
ceed in the following way. We assume that the full ampli- 
tude M12 is the sum of the NP amplitude A NP e^ NP and 
of a SM real amplitude containing both short- and long- 
distance contributions, Asm- We take A^m to be flatly 
distributed in the range [—0.02, 0.02] ps -1 , so that it can 
saturate the experimental bound in Tab. [Til an d derive 



Parameter 


Value 


Ref. 




(-0.24 ± 0.43 ± 0.30) • 10" 3 


m 




(-0.20 ±0.41 ± 0.29) • 10~ 3 




y+ 


(9.8 ±6.4 ±4.5) ■ 10~ 3 




/ 

V- 


(9.6 ±6.1 ±4.3) ■ 10~ 3 


m 


X 


(7.9 ±3.4) ■ 10~ 3 


[10,i2] 


y 


(3.4 ±2.8) ■ 10~ 3 


[10,i2] 


0(°) 


-12 ± 18 


[10,^2] 


I«/Pl 


0.88 ±0.31 


[10] 


yep 


(11.2 ± 3.2) ■ 10~ 3 


[9, 13} 




(-1.7 ±3.0) ■ 10~ 3 


[9, 13] 


COS Sk-k 


1.09 ±0.66 


□jj 


Rm 


(21 ± 11) ■ 10 -3 


\U] 


T D (ps) 


0.4101 ± 0.0015 


[16] 



TABLE I: Experimental results used in our analysis. For 
^4r, yep and Rm we have used the Heavy Flavor Averaging 
Group (HFAG) averages as of May 2007. For x, y and <f> we 
have performed our own combination of experimental results 
as the HFAG averages are obtained assuming no CP violation. 



Parameter 


68% prob. 


95% prob. 


X 

y 

Sktt 
If 1 - 1 

|Afia| (ps" 1 ) 

$12 (°) 
|F 12 (ps- 1 ) 


(6.2 ±2.0) • 10 -3 
(5.5 ±1.4) ■ 10- 3 
(-31 ±39)° 

(1±7)° 
-0.02 ± 0.11 
(7.7 ±2.4) ■ 10~ 3 
(2 ± 14) U (179 ± 14) [ 
(13.6 ±3.5) ■ 10" 3 


[0.0022, 0.0105] 

[0.0027, 0.0084] 
[-103°, 28°] 
[-15°, 17°] 
[-0.27,0.25] 

[0.0030, 0.0127] 
-30,36] U [144,210] 

[0.0068, 0.0207] 


TABLE II: Results on mixing and CP 


violation parameters. 




68% prob. 


95% prob. 


A N p (ps -1 ) 
<?W (°) 

[- 


[0, 0.006] [0, 0.02] 
[-180, -149] U [-180, -112] U [-68, 180] 
-31,39] U [141,180] 



TABLE HI: Allowed ranges for the NP amplitude. 



from the $i 2 vs |Mi 2 | distribution the p.d.f. for Anp 
vs </>np, barring accidental order-of-magnitude cancella- 
tions between SM and NP contributions. The results, 
reported in Tab. IIIII and shown in Fig. [21 provide a new 
constraint that should be fulfilled by any extension of the 
SM. We see that the lack of knowledge of the SM contri- 
bution causes a dilution of the bound on </>np ■ Clearly, if 
a reliable estimate of Asm were available, the constraint 
would be much more effective. Notice also that if |Mi 2 | 
is dominated by NP, then <^np ~ < &i2 and the NP phase 
can be experimentally accessed. 
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FIG. 1: Probability density functions of the combined fit from 
Tab. HI Dark (light) regions correspond to 68% (95%) proba- 
bility. 
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FIG. 2: Probability density functions of the combined fit from 
Tab. [J projected ont/vsi (top left), <f> vs \q/p\ — 1 (top right), 
$12 vs IM12] (center left), Anp vs <^np (center right), Anp 
(bottom left) and <^>np (bottom right). Dark (light) regions 
correspond to 68% (95%) probability. 
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350 


350 


0.033 


0.0056 


0.0020 


500 
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0.049 


0.0080 


0.0029 


1000 


1000 


0.10 


0.016 


0.0062 


500 


1000 


0.14 


0.011 


0.0044 


500 


350 


0.032 


0.0068 


0.0025 



Si = 0.87 ± 0.03 B 2 = 0.82 ± 0.03 B 3 
Bi = 1.08 ± 0.03 B 5 = 1.46 ± 0.09 



1.07 ±0.09 



TABLE V: Upper bounds at 95% probability for |(<5f2) AS | 
for various values of squark and gluino masses (in GeV). 



TABLE IV: B parameters defined as in ref. [21( interpolated at 
the physical D meson mass, renormalized at the scale /1 = 2.8 
GeV in the Landau- RI scheme. 



We now turn to the MSSM and consider the bounds on 
(^12)45 that can be obtained from the determination of 
Anp and 0np discussed above. To this aim, we focus on 
gluino exchange and use the full Next-to-Leading expres- 
sion for the Wilson coefficients [19( and for the renormal- 



ization group evolution down to the hadronic scale of 2.8 
GeV [20|. For the matrix elements, we extrapolate the 
results of ref. 21 1 as given in Table IIVI (see also ref. [22| 
for another recent calculation of Bi). 

To select the allowed regions on the R.e(<5" 2 ) j4s - 
Im((j>" 9 ) A p planes, we use the method described in 
ref. [23|. The results are presented in Fig. [3] for a refer- 
ence value of 350 GeV for squark and gluino masses. We 
consider three cases. First, a dominant LL mass inser- 
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tion. The case of a dominant RR insertion is completely 
identical. Second, a dominant LR insertion. In this 
case, chirality-flipping four-fermion operators are gen- 
erated. These operators are str ong ly enhanced by the 
renormalization group evolution 24] , so that these mass 
insertions are more strongly constrained than LL or RR 
ones. Constraints on RL insertions are identical. Finally, 
we can switch on simultaneously {S^ll = (^12) rr- ^ n 
this case, we also generate chirality-flipping operators, 
so that the constraint is much stronger than the case in 
which (6? 2 ) LL » (Jj» 2 ) 

RR- 

In Table fVl we report the bounds on the absolute value 
of the mass insertions for several values of gluino and 
squark masses. Our bounds are typically a factor of ~ 3 
more stringent than those of ref. [6| . 

It is very interesting that SUSY models with quark- 
squark alignment generically predict (Si 2 ) LL ~ 0.2 0. 
We conclude that, to be phenomenologically viable, they 
need squark and gluino masses to be above ~ 2 TeV. 
Therefore, they probably lie beyond the reach of the 
LHC. 

In this Letter, we have analyzed the first experimen- 
tal evidence of D-D mixing recently obtained by the 
BaBar and Belle collaborations. Combining the exper- 
imental results we obtained new constraints on the mix- 
ing amplitude and on NP contributions. We have then 
considered the MSSM and derived new bounds on off- 
diagonal squark mass terms connecting up and charm 
squarks. Finally, we have briefly commented on the 
impact of these new constraints on SUSY models with 
quark-squark alignment. 
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